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ABSTRACT 

We studied the effect of interstellar gas conditions on global galaxy simulations by considering three 
different models for the ISM. Our first model included only radiative cooling down to 300 K, our second 
model added an additional background heating term due to photoelectric heating, and our third model 
uses an isothermal equation of state with a temperature of 10 4 K and no explicit heating or cooling. 
Two common prescriptions for star formation are implemented in each case. The first is based on 
cosmological simulations with a low threshold for star formation but also a low efficiency. The second 
assumes stars form only in high density regions but with a higher efficiency. We also explore the 
effects of including feedback from type II supernovae. We find that the different ISM types produce 
marked differences in the structure of the disk and temperature phases present in the gas, although 
inclusion of feedback largely dominates these effects. In particular, size of the star-forming clumps 
was increased both by background heating and by enforcing an isothermal ISM. We also looked at the 
one dimensional profiles and found that a lognormal PDF provides a good fit for all our simulations 
over several orders of magnitude in density. Overall, despite noticeable structural differences, the star 
formation properties in the disk are largely insensitive to ISM type and agree reasonably well with 
observations. 

Subject headings: galaxies: spiral, galaxies: ISM, galaxies: evolution, methods: numerical, ISM: 
structure 



1. INTRODUCTION 

Star formation is one of the most perplexing processes 
in the galaxy due to its immense complexity on the very 
small scales and its apparent simplicity on the very large. 
What is difficult to know is to what extent the small 
scale physics can be ignored when considering the global 
evolution of galaxies. 

Observationally, we know there is a simple relationship 
between the surface density of gas in a disk galaxy and 
the surface density of star formation (|Kennicuttj )l998l 
ll989HSchmidtlll959h . This should imply that star forma- 
tion is a straight forward, universal process depending on 
average properties of the gas over many kiloparsecs. The 
(perhaps naive) picture is one of gravity acting to collapse 
the gas into giant molecular clouds out of which stars 
form, a process which is hindered by rotational shear, 
thermal pressure, turbulence, magnetic fields, c osmic ray 
press ure and energy injected from supernovae (jMacLowl 
l200l . Yet when we start looking at star formation on 
small, parsec, scales, we find the gas out of which stars 
form to be a turbulent, multiphase medium, strongly af- 
fected by local conditions and not at all indicative of a 
global law. 

So complex is this interstellar medium (ISM) that sim- 
ulations which model it have been forced to consider only 
small section s of the galaxy t o achieve the required res- 
olution (e.g. ISlvz et all 120051: Ucmng fc Mac Lowl 120051: 
Ide Avillez fe Breitschwerdtl |2004[ ) . Until recently, sim- 
ulations which model the entire galaxy have been forced 
to simplify the structure of the ISM to an isothermal 
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or fixed two- or three-pha se body for the c a lculation 
to be at all feasible (e.g. Harfst et al. 200 6| ILTet alJ 
2005t iRobertson. Yoshida. Springel. fc Hernqulstl l2004t 
Semelin fc Combesl l2002f )~ Not that such assumptions 
about the structure of the IS M are groundless. Analyti - 
cal calculations performed by iMcKee fc Ostrikerl (|1977l ) 
put together the now traditional picture of a three-phase 
ISM, largely filled by hot gas from supernovae e xplo- 
sions. Others have suggested modifications (e.g. ICoxl 
[20051 : IN orman fc Ikeucha[l989h but this basic scenario is 
still with us. 

Yet while both of these approaches have been highly 
informative about the evolution and structure of galaxies 
they each have disadvantages. The small box simulations 
are unable to model global properties such as star forma- 
tion histories, disk structure and the Schmidt law, result- 
ing in their main comparison points being restricted to 
galaxies where we are able to observe and measure the 
ISM. Global models, on the other hand, have been un- 
able to compute the evolution of the interstellar medium 
and the effect of this simplification on the disk galaxy's 
own evolution is unknown. 

Recently, however, simulations have been developed 
that are bridging this gap. While still short of the resolu- 
tion attained in kiloparsec-sized simulation boxes, these 
models do achieve the refinement needed to model a 
complete galaxy disk that includes a multiphase ISM. 
This pa per continues our work began in Tas ker fc Brvanl 
(|2006al) (hereafter TB06) which examined a three- 
dimensional isolated galactic disk using an adaptive- 
mesh refinement technique. In that paper, we found that 
a multiphase medium with a large variety of tempera- 
tures and densities, in good agr eement with simu l ations 

-rnanl (l200lD 
Wada et al.l 
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(2002). These results supported the small box simula- 
tions in implying that the traditional view of the inter- 
stellar medium as a three-phase structure in strict pres- 
sure equilibrium has some validity but is a significant 
over-simplification. 

What is less clear is how important this is. The ex- 
istence of the Schmidt law might imply that the de- 
tails of the ISM are not necessary to achieve accurate 
modeling in global galaxy simulations. And indeed, the 
Schmidt law and the observed star fo rmation thresh- 
old in disk galaxies (|Kennicuttl 11989ft has been suc- 
cessfully reproduced both in models with a multiphase 
ISM (TB06) and those with a fixed isothermal ISM 
(|Li, Mac Low fc Klessenll2005[) . On the other hand, the 
observed global structure of the ISM is far from an unin- 
terrupted pool of gas. In addition to its turbulent nature, 
observation of the ISM in the Large Magellanic Clouds 
show a com plex series of HI filaments riddled with holes 
and shells (|Kim et al.l fl998). Two of our neighboring 
spiral galaxies, M31 and M33 show nume rous holes 40 pc 
to 1 kpc wide (|Deul fc den Hartod [1990) , and our own 
Milky Way produces plum es of gas that rise off the disk's 
surface (|Otte et al.ll2003h . Not only can this not be mod- 
eled without global multiphase ISMs, it seems impossi- 
ble that this structure cannot play an integral part in 
the disk's evolution. Moreover, the role of stellar feed- 
back remains something of a mystery in both isolated 
galaxy simulations and in cosmological runs where in- 
correct modeling is frequently cited as a possible cause 
of disc repancy with observational data (jTasker fc Bryanl 
2006b, and references therein). Improved modeling of 
the ISM on a global scale could result in a fuller under- 
standing of feedback, ultimately allowing more accurate 
feedback routines in large cosmological simulations where 
resolution of the individual galaxies is not yet possible. 

In this paper we will compare global models of isolated 
disk galaxies with three distinct ISM types. The first of 
these will include radiative cooling down to 300 K, the 
second type will contain radiative cooling and a back- 
ground photoelectric heating source and our final type 
will have an isothermal ISM at a constant temperature 
of 10 4 K. For each of these models, we will test two com- 
mon prescriptions for star formation and the effects of 
feedback from type II supernovae. The resulting struc- 
tures are contrasted and the star formation properties of 
the disks compared with observations. 

For these simulations, we use a high-resolution adap- 
tive mesh refinement (AMR) code which includes a full 
treatment of self-gravity of the gas rather than the fixed 
potential which is often used and a more sophisticated 
treatment of star formation and feedback. Our simu- 
lations concentrate on hydrodynamical effects, ignoring 
magnetic fields and cosmic ray pressure. 

In section 2, we describe our computational approach, 
including details of the code we are using and lay out 
the initial conditions for the problem. Sections 3 and 
4 will focus on the disk structure and the properties of 
the ISM while section 5 will look at comparisons with 
observations. 

2. NUMERICAL METHODS 

2.1. The Code 

The simulations were performed using the hydrody- 
namics adaptive mesh refinement (AMR) code, Enzo, 



described in iBrvan fc Normanl <l!997f): iBrvanl (l!999f) 



iNorman fc Brvanl (ll999h ; lBrvan. Abel fc Normail (|2001l ) 
and lO'Shea et al I (|200l . The AMR technique is par- 
ticularly strong in this work where the resolution of 
a complex multiphase medium is of paramount impor- 
tance. The grid cells form natural boundaries which al- 
lows g as with a range of temperatures to coe xist and 
evolve (|Slvz et al.ll2005l : iTasker fc Bryanl l2006al ). Other 
codes using particle-based techniques frequently have an 
over-mixing problem, causing unp hysical radiative losses 
unless algorithmic steps are taken (|Marri fc Whitdl2003l: 
ISpringel fc Hernqui st 2003). As a result, the majority of 
previous simulations have been unable to properly model 
a multiphase interstellar medium, making them unable 
to access its importance in star formation. 

For these simulations, we use a three-dimensional peri- 
odic box of side 1 h -1 Mpc. With a root grid of 128 3 and 
8 levels of refinement; our smallest cell size, and there- 
fore maximum resolution, is approximately 50 pc. For 
our higher resolution run, this was decreased further to 
25 pc. 

To evolve the gas through time, Enzo used a three- 
dimens ional version of t he ZEUS hydrodynamics algo- 
rithm (| Stone fc Normanl 1 1992f ) Radiative gas c ooling 
followed the cooling curve of iSarazin fc White! (|l987fr 
down to temperatures o f 10 4 K and then rates from 
iRosen fc Bregmanl (|1995f ) down to 300 K. (The excep- 
tion to this are for the simulations performed with an 
isothermal equation of state for the gas, where no cool- 
ing is allowed). This bottom temperature threshold is 
still above what would be foun d in th e dense molecular 
clouds, but IRosen fc Bregmanl (|1995f) argued that this 
crudely compensates for physical processes not modeled, 
such as magnetic fields, turbulence, and cosmic-ray pres- 
sure. This temperature range does take us to th e uppe r 
limit of the cold neutral medium IWolfire et alJ ([2003) , 
allowing us to sample a realistic spread of phases in the 
gas. 

In addition to cooling, the gas can also be heated 
through supernovae feedback (described below) and pho- 
toelectric heating. Photoelectric heating, that is the 
photo emission of UV-irradiated dust grains, is thought 
to be the dominant factor i n the formation of the cold 
and warm neutral mediums (jWolfire et al. 1995). If pho- 
toelectric heating is turned on, Enzo includes the term 
Tpe = 5.1 x 10~ 26 ergs _1 to the energy equation for the 
gas in a scheme based on Uoung fc Mac Lowl (|2005f ) , but 
without the dependence on the height above the disk 
adopted in that work. The value of Tpe is dependent on 
the incident radiation field and, as such, could be tied 
with the star formation rate in the gas. However, since 
we are considering Milky Way -sized galaxy disks, we fol- 
low [joung^&^a£jx>w| (|2005[ ) and adopt a number con- 
sistent with the local interstellar value. Since it is likely 
that this value was higher at earlier times, we also per- 
form a run with a heating source of r^g = 1.41 x 10~ 24 
ergs -1 = 30rpe, a value that balances the cooling rate at 
densities of 1.0 cm~ 3 at our initial temperature of 10 4 K. 
However, the increase in heating makes a relatively small 
difference to the star formation properties of the disk, so 
this run is not included in our main analysis but added 
to the discussion only. Note that we do not include pho- 
toionization heating, which would require full radiative 
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transfer to model properly. 

Enzo will form a star parti cle in a grid cell if 
it fulfills the follo wing criteria (|Cen fc Ostrikerl Il992t 
lO'Shea et~aT1 12004) : (i) the baryon density in the grid 
cell exceeds a designated threshold density, (ii) the mass 
of gas in the cell exceeds the local Jeans mass, (iii) there 
is convergent flow into the cell (i.e. V • v < 0) and (iv) 
the cooling time of the gas in the cell is less than its dy- 
namical time (t coo j < ), or the gas is at the mini- 
mum temperature allowed of 300 K. We will also consider 
one set of galaxy disks which are modeled with a purely 
isothermal interstellar medium. In this case, the above 
criteria are relaxed to only include the first three rules, 
since the gas is unable to cool. 

If a grid cell meets the required criteria for star forma- 
tion, gas is removed from the cell and a star particle will 
be formed with mass calculated by: 

At n 

m* = e- pgasAz (1) 

t dyn 

where e is the star formation efficiency (more properly 
the efficiency per dynamical time) , At is the size of the 
time step, t^y n is the time for dynamical collapse and 

Pgas is the gas density. A final, purely computational, 
criteria is that, even if a grid cell fulfills all the above cri- 
terion, a star particle will still not be formed if its mass 
is less than a given minimum value of (for this work) 
10 4 M . Since the reason for this clause is purely numer- 
ical (a large number of small stars would greatly slow 
down the simulation) an override exists that creates a 
star below the minimum mass with a probability equal 
to the ratio between the mass of the would-be star par- 
ticle and the minimum star mass. If this occurs, the 
resulting star mass is that of the minimum star mass or 
80% of the mass in the cell, whichever is smaller. The 
star particle's formation is then spread out over roughly 
a dynamical time, to mimic the formation of stars in the 
giant molecular clouds who follow the same time scale. 

The star particles themselves are modeled as a colli- 
sionless system using the N-body method. They gravita- 
tionally interact with the gas by mapping their positions 
onto the grid via a cloud-in-cell technique to produce a 
discretized density field. The number of star particles 
formed varies greatly depending on the values chosen for 
the threshold density and star formation efficiency de- 
scribed above. For our 'C'-type star formation method 
where these values are low (see section I2.3j) , roughly five 
million particles were used by the end of the simulation. 
In our 'D'-type algorithm, these values are higher and 
the number of star particles formed is reduced to around 
40,000. 

Enzo can also include stellar feedback from type II su- 
pernovae, often suggested as the main driving force for 
self-regulated star formation. When this feedback op- 
tion is used, as it is in about half our simulations, then 
10~ 5 of the rest-mass energy of generated stars is added 
to the gas' thermal energy over a time period equal to 
^dyn' This is equivalent to a supernova of 10 51 erg for ev- 
ery 55 M Q of stars formed. The energy is deposited into 
the gas over one dynamical time or 10 Myrs, whicheve is 
longer. During this period the energy is injected into the 
cell closest to the particle's current location. 

For all simulations, we adopt a cosmological model of 



TABLE 1 

Overview of simulations performed 







min iVr 
(vc) 


6 


^thresh 

( cm -3 ) 


r 


Iso 






f c, 


50 


0.05 


0.02 


No 


No 


No 


ISM #1 J 


I D 


50 


5 


10 3 


No 


No 


No 




| CFDBCK 


50 


0.05 


0.02 


Yes 


No 


No 




[ DFDBCK 


50 


0.5 


10 3 


Yes 


No 


No 




f HC 


50 


0.05 


0.02 


No 


No 


Yes 


ISM #2 J 


HD 


50 


0.5 


10 3 


No 


No 


Yes 




HCFDBCK 


50 


0.05 


0.02 


Yes 


No 


Yes 




HDFDBCK 


50 


0.5 


10 3 


Yes 


No 


Yes 




I HDHIRES 


25 


0.5 


10 3 


No 


No 


Yes 




r 10 


50 


0.05 


0.02 


No 


Yes 


No 


ISM #3 | 


ID 


50 


0.5 


10 3 


No 


Yes 


No 




{ ID JEANS 


50 


0.5 


10 3 


No 


Yes 


No 


a ACDM universe with f2 


m = 0.3, 


= 0.7 and H 







67kms~ 1 Mpc -1 . 

2.2. The Initial Conditions 

The initial conditions for our disk are the same as the 
simulations performed in TB06, with an ideal isothermal 
gas disk of temperature of 10 4 K and 7 = 1.67 whose 
density profile is given by 

p(r,z) = p e-^sech 2 ^^, (2) 

sitting in a static dark matter halo. The major difference 
between these initial conditions and the ones set out in 
TB06 is the gas mass in the disk which we choose to 
be 6 x 10 10 M Q , six times higher than in TB06, bring- 
ing it into line with the estimated total disk mass of 
the Milky Way (we do not begin with any stars). We 
also change the disk dimensions slightly from TB06, se- 
lecting a scale radius r = 3.5 kpc and a scale height 
zq = 325 pc. Together with the gas mass, these choices 
fix the value for po = O.6M pc -3 . Dark matter is in- 
cl uded as a static halo potential in the form described 
bv lNavarro. Frenk fc White! (| 19971) . The disk is initially 
borderline stable, with a Toomre Q parameter (described 
in section 15. 2p of 0.5 at the center of the disk rising to 
10 at the edge, but quickly cools to fragment. This is 
discussed in more detail in TB06. 

2.3. Summary of Performed Runs 

Our simulated galaxy disks are divided into three cat- 
egories depending on the nature of their ISMs. For the 
first four simulations listed in table [TJ the gas is allowed 
to radiatively cool down to 300 K via the cooling curves 
described in section 12.11 This is the same set-up (al- 
though for a heavier disk) that we presented in TB06. 
The next five runs also allow radiative cooling along the 
same curve, but include an additional photoelectric heat- 
ing term, as described in section [2~T1 The last three sim- 
ulations use an isothermal equation of state, a popular 
assumption in global disk models where it has been hard 
to resolve a multiphase medium. The temperature of all 
the gas in these runs is fixed at 10 4 K. 

For each of these three ISM models, we consider two 
different star formation routines. These were presented 
in TB06 and we keep the same notation of 'O and 'D' 
type as used in that paper. C-type star formation has 
an efficiency appropriate to the galactic disk as a whole, 
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and one typically used in large-scale cosmological simu- 
lations (5%). Because this is a global average, we use a 
low density threshold, allowing stars to form in relatively 
low density regions (providing they meet the criteria out- 
lined in section 12. ip . Since this means we do not follow 
the formation of the densest clumps, we use a Schmidt- 
like law to model the star formation rate. D-type star 
formation, on the other hand, makes use of the resolu- 
tion in the disk and confines star formation to the dens- 
est structures which, at our resolution, compare to the 
largest giant molecular clouds. Stars therefore form only 
where the density is high (10 3 cm~ 3 ) but with a much 
higher efficiency per dynamical time of 50%. 

Simulations which include stellar feedback from type 
II supernovae are performed for both the disks with ra- 
diative cooling (CFDBCK and DFDBCK for C-type and 
D-type star formation respectively) and for disks with 
radiative cooling and background heating (HCFDBCK 
and HDFDBCK). 

In addition to these simulations, we also perform runs 
designed to test the robustness of our results. HDHIRES 
is the same as HD, but with a root grid which has twice 
the spatial resolution and so eight times the mass res- 
olution. ID JEANS uses identical conditions to ID, but 
adds another refinement criteria, resolving a cell when 
the Jeans length drops below four cell widths (at least 
until we r each the maximum re finement level), as sug- 
gested by iTruelove et all (JT997) . As mentioned in sec- 
tion 12.11 we also perform a run with a higher heating 
rate. This is discussed in the discussion section at the 
end of this paper. 

The run times for these simulations depended greatly 
on whether stellar feedback was included. Without feed- 
back, a run typically took ~ 40 hours on 8 processors of 
an Opteron Beowulf cluster. When feedback was intro- 
duced, the run took ~ 130 hours on 64 processors of a 
Xeon cluster. 

3. THE STRUCTURE OF THE DISK 

During the initial evolution of our galaxy disks, we see 
the gas fragment, triggering a starburst. This process is 
described in detail in TB06, but for this work we want 
to avoid this primary collapse and compare our disks 
when they have reached a settled state but before any gas 
depletion becomes an issue. To judge where this point 
is, we examine the gas density in the disk over time. 
In the first 200Myrs, we see a clear and rapid drop in 
density corresponding to the initial fragmentation. The 
density then settles down, decreasing slowly as stars are 
formed until after around 500 Myrs were the profile starts 
to deteriorate, especially in the non-feedback cases, as 
the majority of the gas is now converted to stars. Our 
images and plots are therefore largely taken between 200 
- 400 Myrs when the initial starburst is over and where 
the disks still have a substantial supply of gas. 

3.1. Imaging the ISM 

The global structure of the ISM can be viewed using 
face-on projections of the disk as shown in Figures [IJ 

and [3] Figure [1] shows projections for our first ISM 
type, the top four simulations listed in table [TJ all of 
which include radiative cooling but no background heat- 
ing. From left to right, the images show gas density, 
temperature, pressure and stellar density. The top two 



columns which shows runs without feedback clearly show 
the region where the gas is gravitationally unstable. As 
in TB06, we see that the central gas has collapsed to form 
dense filaments and knots out to a clearly defined radius. 
Outside this region, the gas is stable and remains unper- 
turbed. Earlier images show the same pattern forms for 
the initial evolution of all the runs, with a circular per- 
turbation moving outwards from the center of the disk 
which then collapses tangentially to form filaments and 
knots of gas. The introduction of feedback (bottom two 
columns) destroys these filaments, smoothing the gas dis- 
tribution regardless of the type of star formation used. 
The effect of this smoothing is shown both in the tem- 
perature distribution and with the disk's pressure. Hot 
outflows are seen and cold gas is no longer confined to 
the dense star forming knots, but rather is blown around 
the disk. We will see later than this acts not only to 
change the balance of the phases in the ISM but also to 
suppress the star formation as knots of gas are destroyed 
before they can collapse into stars. It also affects the 
stability in the outer regions of the disk which we will 
discuss more thoroughly in section f5. 21 

The feedback temperature projections also show us the 
first evidence of the different effects of the C- and D-type 
star formation routines. The C-type stars cause rela- 
tively even outflows (something we will return to in the 
next section), resulting in a more uniform disk. The D- 
type, on the other hand, concentrates the outflows, and 
so we see one side of the disk is much colder than the 
other side. The position of these outflows changes as the 
disk evolves, but the separation of the phases is always 
present. This result is unsurprising if we consider the 
physics of each stellar type; the D-type stars only form 
in the densest regions, whereupon they form efficiently 
and in large numbers. This produces a focusing of stellar 
material which is transferred to a concentrated injection 
of energy from supernovae. On the other hand, C-type 
stars form more uniformly, producing a smoother distri- 
bution and hence an even injection of energy. This is 
shown most clearly in the final column displaying stellar 
distributions. The C-type stars, forming at much lower 
densities, extend smoothly out until the threshold ra- 
dius for gravitational collapse. The D-stars, by contrast, 
form only in the densest areas of the disk, confining them 
to the central region where the gas has collapsed both 
radially and tangentially into dense knots. This result 
differs somewhat from what was found in TB06, where 
the lighter disk meant that both C- and D-type star for- 
mation algorithms produced stars only in the completely 
collapsed regions, whereas here we see C-type stars form- 
ing in the mildly perturbed areas of the disk (see also 
section I5T2")) . 

The feedback outflows also disturb the pressure in the 
disk. Non-feedback runs C and D show small-scale vari- 
ations in the pressure distribution, with the dense, star 
forming knots of gas being at a higher pressure than the 
surrounding ISM. The inclusion of feedback, especially in 
the focused D-type star formation, upsets this, showing 
the large-scale hot outflows to be over-pressurized with 
respect to the disk. 

Our second ISM type includes background heating in 
addition to radiative cooling. Figure [2] shows the pro- 
jections for runs HC, HD, HCFDBCK and HDFDBCK 
for the same properties as Figure [T] Unlike the first ISM 
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mS 14.2 11.8 9.50 2.S4 4.32 5.73 7.27 -15.4 -1 3.0 -10.6 -a .26 9.00 11.9 14.S 17.7 

Fig. 1. — Projections at 377Myrs of (left to right) gas density, temperature, pressure and stellar density for the runs with ISM #1 which 
include cooling but no photoelectric heating. Top to bottom, simulations shown are C, D, CFDBCK, DFDBCK. Images are 60kpc across. 
All scales are to the base-10 logarithm, and gas and star particle density is measured in MqM/dc - 2 , temperatures in K, and pressure on 
an arbitrary scale. 



type, we see a notable difference in the C and D star for- 
mation types in the runs without feedback. The C-type 
star formation produces large voids that contain hot, low 
density gas. This shape comes out less clearly in the 
stellar distribution, but the dense star clusters are con- 
fined to the filaments surrounding these holes. The holes 
have largely vanished in the disks that include feedback, 
although they all display a porous nature. The same 
holes are also seen in run C, but only during the first 



250Myrs of the simulation. After that point, the circu- 
lar wave speeds up with respect to the HC simulation, 
and decreases in strength. Exactly why this occurs is 
not clear. One explanation is that the extra heating in 
the HC case provides an added pressure that stabilized 
the circular mode. Alternatively, this is a numerical ef- 
fect and with improved resolution the porous structure 
would be retained in the C case as well as the HC sim- 
ulation. Either way, the presence of holes in the disk is 
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Fig. 2. — Projections at 377Myrs of (left to right) gas density, temperature, pressure and stellar density for the runs with cooling and 
photoelectric heating. Top to bottom, simulations shown are HC, HD, HCFDBCK, HDFDBCK. Images are 60kpc across. All scales are 
to the base-10 logarithm, and gas and star particle density is measured in M c(otMpc~ 2 , temperatures in K, and pressure on an arbitrary 
scale. 



an interesting event and one that has been seen both in 
our own galaxy and, more dramatically, in the HI map 
of the LMC. Their presence has traditionally been put 
down to stellar winds and supernovae exp losions evac- 
uating the cool IS M (jvan der Hulstl 1 1996T ) . However, 
iRhode et all (flflfll was unable to find evidence of rem- 
nant star clusters in the center of the HI holes in the ir- 
regular galaxy Holmberg II, suggesting that supernovae 
were not present there. The issue was also investigated 



th eoretically via two- dimensional simulations performed 
bv lWada et alj (|2000D of an LMC-type galaxy which sug- 
gested that gravitational and thermal instability alone 
are enough to create a porous ISM and that these are ac- 
tually disrupted in the presence of frequent supernovae. 
This would appear to agree well with our findings that 
the non- feedback run contained holes, but these largely 
vanish when feedback is included in the simulation. 
The run HD with D-type star formation also shows 
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Fig. 3. — Projections at 377Myrs of (left to right) gas density, pressure and stellar density for the runs with an isothermal ISM. 
Temperature is 10 4 K throughout simulation box. Simulations shown are IC (top) and ID. Images are 60kpc across. All scales are to the 
base-10 logarithm, and gas and star particle density is measured in M doiMpc~ 2 , temperatures in K, and pressure on an arbitrary scale. 



some evidence of cavities, but not as strongly as in HO 
What is noticeable is that the dense knots of gas are far 
more evident in this image than in the equivalent run D 
for ISM #1 in Figure [T] The overall filament structure in 
both HC and HD is reduced compared to the non-heated 
simulations, suggesting that the effect of heating is to act 
against the collapse, increasing the Jeans length to allow 
only the larger perturbation to form knots. These dense 
knots, however, extend out further in the heated case 
than the non-heated case, but the smaller perturbations 
are smoothed out. We will return to this quantitatively 
when we consider disk stability in section [5T2l 

Feedback again acts to smooth the gas distribution, al- 
though less effectively than in the CFDBCK and DFD- 
BCK cases, allowing dense knots of gas to survive, again 
indicative of them being larger and more tightly bound 
than in the ISM #1 case. We see the same symmetry to 
the feedback outflows with C-type star formation as with 
Figure [T] and the focused ejections of the D-type simu- 
lations, where the star formation is confined to a much 
smaller area. 

The pressure distribution in the disk is almost entirely 
isobaric for the runs without feedback, with only the 
densest knots in the D-type star formation in HD be- 
ing at a slightly higher pressure. The heating of the low 
density background gas also raises its pressure and causes 
it to be in pressure equilibrium with the disk. Feedback 
again acts to disrupt this in the disk, most markedly 
in the D-type star formation case where the hot gas is 
marked by considerably higher pressures. 

Our final ISM type is for the disks with an isothermal 
ISM at a constant temperature of 10 4 K. The projections 



of these runs are shown in Figure [3J We see the effects 
of stopping the disk from cooling in the size of the knots 
of gas, which are visibly larger than in either ISM #1 
or ISM #2 due to the higher Jeans length throughout 
the whole disk. As a result, the fragmentation is limited 
to the disk's central region which, especially in the case 
of the C-type star formation, becomes depleted of gas. 
The resultant cavity looks similar to that found in HC 
although it is smaller and corresponds to regions of dense 
star clusters, rather than a lower density of stars. In this 
respect, the C-type star formation looks more like the D- 
type, with the majority of star formation in large clusters 
in the disk center. Later evolution shows this central 
cavity growing, leaving a gas deficit void. Overall, the 
disk takes on a smoother appearance. 

The production of larger, if fewer, dense gas knots has 
a dramatic effect on the stellar density for the C-type star 
formation run. A nalysis of the star clusters (described in 
iGill et all (|2007| )) produced by the large knots of gas re- 
veal masses up to I0 10 M with 1/2 mass radii of order a 
cell size. These huge dense clumps gravitationally inter- 
act with close star particles, accelerating a small fraction 
of them up to 1000 kms -1 resulting in ejection from the 
disk in high velocity streams. By contrast, the largest 
star clusters formed in the C run are a factor of ten less 
in mass and much more diffuse resulting in a significantly 
smaller gravitational pull on nearby star particles. In the 
case of D-type star formation, the star clusters are more 
tightly bound (since they form only in the densest gas) 
which makes it harder for star particles to be ejected at 
high velocities. 

Outflows of this magnitude and speed in the IC case are 
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Fig. 4. — Edge-on projections of the baryon density in the 
CFDBCK and DFDBCK simulations after 142 Myrs. Images are 
~ 210 kpc across. Both simulations include feedback and radia- 
tive cooling, but the left-hand image (CFDBCK) has a low density 
cut-off and low efficiency for star formation whereas the right-hand 
image shows the disk with a high density cut-off and high star for- 
mation efficiency. 
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Fig. 5. — Average mass flux out of the disk between t = 283 - 
472 Myrs for runs CFDBCK and DFDBCK. Solid line shows the 
outflow from the disk only while the dashed line shows the net flow 
above and below the disk. 

not observed in real galaxies, and this is a point against 
using the isothermal ISM model with the C-type star 
formation algorithm. However, the ejection of star par- 
ticles by cluster interactions raises an interesting ques- 
tion about star formation in the outer parts of the disk. 
Observational ly, low luminosity stars have been foun d 
at large radii ([Ferguson et alJl2002t iBoissier et alJ l2006) , 
beyond the point where the disk is traditionally gravi- 
tationally stable. How these stars got there is an open 
question but they may have been produced during satel- 
lite interactions, or they may have been thrown there 
from interactions within the disk. Figure [3] suggests that 
interactions between heavy star clusters could potentially 
produce this effect. 

The pressure projections show the isothermal model to 
be the least isobaric of the three ISM types. This is un- 
surprising when we consider that fixing the temperature 
forces the pressure to mirror the density distribution. We 
see here that the gravitationally collapsed structures are 
over pressurized with respect to the disk and voids of gas 
have low pressure. 

3.2. The Vertical Profile 

The vertical structure of the disk and — linked with 
this — the interaction between the disk and the halo, 
are tho ught to play a vital part in the di sk's evolution. 
Indeed, Ide Avillez fc Breitschwerdtl (|2004f ) cite the cycle 
of gas between the disk and the halo as being one of the 
primary factors in determining the phases of the ISM 
in their small-box simulations. They find the ability of 



hot gas to lift off the disk's surface acts as a pressure 
release valve on the ISM. This i s in k eeping with the 
work done bv lNorman fe Ikeuchil (|1989f ) whose model of 
the hot gas being removed vertically through galactic 
chimneys all owed a lower filling f a ctor t han in the original 
picture from iMcKee fe Qstrikerl (|1977Tj , more in keeping 
with observations. 

Figure [4] shows the vertical projections of the bary- 
onic density in our feedback runs for ISM #1, CFD- 
BCK and DFDBCK (the results are similar for HCFD- 
BCK and HDFDBCK in ISM #2). The left-hand im- 
age shows the simulation with C-type star formation, 
whereas the right-hand one has D-type stars. Both galax- 
ies are clearly injecting a significant amount of material 
out of the disk and into the halo, consistent with both 
iNorman fe Ikeuchl (Il989f) analy t ical p icture of the ISM 
and Ide Avillez fc Breitschwerdtl (|2004f ) small box simu- 
lations. We have already seen in Figures [1] and [2] that 
the outflow gas is highly pressurized, in keeping with the 
release valve idea. This is in contrast with the simu- 
lations without feedback, where the gas is confined to 
the disk's surface. These images are taken shortly after 
the outflows begin in the disk, at 142 Myrs where the 
difference between the C- and D-type star formation is 
sharply apparent: D-type star formation not only results 
in asymmetrical distribution of gas within the disk as 
seen in section I3TTI but the gas is also ejected unevenly 
from the disk's surface. This is most likely a result of the 
density-threshold in the D-type star formation prescrip- 
tion, which can only be met in large clumps and so is 
naturally highly inhomogeneous. The right hand panel 
of Figure [4] shows one outflow extending away from the 
bottom of the disk which has almost reached the image's 
edge. A second outflow, either just starting up or con- 
densing to fall back in, is shown at the top of the disk. 
The C-type star formation by contrast produces an even 
distribution either side of the disk, suggesting gas is be- 
ing emitted from both sides largely at the same time. 

Despite the symmetrical differences in the gas ejec- 
tion of the two star formation types, neither run sees 
the majority of the gas leaving the gravitational pull of 
the galaxy. Rather, the gas cools above the disk and falls 
back down in a galactic fountain effect. This can be seen 
in Figure which shows both the net mass flux as a 
function of height from the disk and also the outward- 
bound only gas. The mass flux is averaged over almost 
200 Myrs, from t = 283-472 Myrs and over this time range 
the overall outflow from both C- and D-type star forma- 
tion is roughly symmetrical above and below the disk. 
The net mass flux, however, is much smaller than the 
outflow, showing that the majority of the gas returns to 
the disk. The size of the outflow is strongly dependent on 
the star formation type. D-type star formation, with its 
concentrated clumps, produces a significantly stronger 
outflow than the C-type, reaching maximum mass fluxes 
of 30M Q yr~ 1 compared to around lOMoyr -1 . In both 
cases, the outflows are largely restricted to only a few 
kpc away from the disk. In the C-type case, there is 
almost nothing outside this region whereas the D-type 
shows evidence of outflows extending further away from 
the disk and infalling gas up to heights of 20 kpc. 

A more quantitative way of studying the vertical dis- 
tribution in the disk is to examine the one-dimensional 
profiles. Figure [6] shows the variation of baryon density 
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Fig. 6. — The vertical density profile at 472 Myrs into the evolu- 
tion of the simulations listed. This is the density profile averaged 
within two radial scale-heights. 
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Fig. 7. — The vertical pressure profile at 472 Myrs into the evo- 
lution of the simulations listed. The profile is also averaged within 
two radial scale-heights. 



with height after 472 Myrs for each of the ISM types. 
The effect of changing the properties of the ISM is very 
striking in this plot, especially for the runs that do not 
include feedback. In the top panel, simulations C and D 
with ISM #1 are shown which have a peaked profile with 
vertical scale height of about 50 pc (and which is prob- 
ably unresolved in the simulation). When background 
heating is added for HC and HD in ISM #2 (middle 
plot), the disk broadens out, producing a scale height of 
roughly twice that in ISM #1. It is worth noting how- 



ever, that our background heating term is a constant. A 
more peaked profile might have been achieved by vary- 
ing the heating rate wit h disk height as was done in 
iJoung fc Mac Low! (|2005l ). The background heating also 
affects the two star formation types in different ways; 
the HD profile is broadened to a greater extent than HC, 
which becomes flattened in the central most region of 
the disk. In the isothermal case, the disk is more insen- 
sitive to star formation type, although the D-type star 
formation in ID extends to a higher z than IC. The pro- 
file shape and scale height are very similar to HD with 
z ~ 100 pc. 

The addition of feedback has the same effect on both 
ISM #1 and ISM #2, producing a far flatter and broader 
profile, indicative of the destruction of the central dense 
gas knots and material being distributed further out from 
the disk by the outflows. Like ISM type, star formation 
type has little impact in this case. 

Our higher resolution simulation, HDHIRES, closely 
follows its lower resolution counterpart, HD. Likewise, 
the run ID JEANS, specifically designed to resolve where 
the Truelove criteria for resolution is not met, is almost 
identical to ID. This is true for all other plots in this 
paper. 

Over time, the vertical profile remains fairly stable, al- 
though runs with C-type star formation decrease in den- 
sity more quickly than the D-type as the gas is used up 
in the disk. Feedback acts to reduce this effect and in- 
deed, it suppresses the star formation as we shall see in 
section [5] Heating the disk also suppresses star forma- 
tion in the HD run, but has markedly less effect on HC 
which continues to lose gas. The exception to this are the 
isothermal runs which are particularly stable, especially 
ID which shows much less evolution over the course of 
the simulation. 

The vertical pressure profile of the disks are shown 
in Figure [7] Most noticeably, runs with ISM #1 and 
ISM #2 have largely isobaric profiles as indicated in the 
projections. The isothermal runs (ISM #3) meanwhile 
have pressure proportional to density and show a clearly 
peaked distribution. The introduction of feedback in 
ISM #1 and #2 raises the pressure, corresponding to 
the high pressure regions we saw in the radial projec- 
tions in Figures Q] and [2 What is apparent now, how- 
ever, is that the disk is in greater pressure equilibrium 
across its height, as gas is ejected from the disk's surface. 
This shows that the gas ejection seen in Figure 0] is in- 
deed acting as a pressure release system, as described in 
iNorman fc Ikeuchil (|1989h . 

4. PHASES IN THE ISM 

4.1. The Density PDF 

Section [3] shows that the different types of ISM mod- 
els are making a noticeable difference to the appearance 
of the galaxy disk. This brings us back to our original 
problem: if the ISM conditions affect the galaxy so much, 
why do observations show such a simple relation between 
local gas surface density and star formation properties on 
lar ge scales? A way to explain this has been suggested 
by lElmegreenl (|2002t ). who noted that if the density dis- 
tribution were everywhere a log-normal distribution, and 
if star formation occurred above a fixed density, then the 
Schmidt law would result. This means that despite its 
complex nature on the very small scales, the fraction of 
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Fig. 8. — Probability distribution function of the volume 
weighted gas density in the simulated galaxy disks at 377 Myrs 
and at 1.32 Gyrs. A lognormal distribution is overlayed in the top 
two panels. 

gas dense enough to form st a rs wo uld always be known 
from this profile. iSlvz et ail (|2005f ) does point out that 
this does not completely hold, since no spatial informa- 
tion would be available from such a PDF, resulting in the 
possibility that the dense regions do not contain enough 
mass to exceed the critical Jeans limit for collapse. That 
aside, the presence of a universal density distribution 
function in our disks could be indicative that such a link 
between the small and large scales is possible. 

Previous simulation work has found evidence that a 
single profile could be fitted to a range of galaxies, but 
struggle to agree on the shape. The debate largely cen- 
ters on whether the PDF is best fitted by a lognormal or 
a power-law cur ve at high densiti es (where star forma- 
tion will occur). iScalo et ail (|1998l ) puts this discrepancy 
down to the nature of the gas, finding that an isother- 
mal flow follows a lognormal distribution, while a gas 
with 7 ^ 1 is better represented by a power-law fit. 
T his view is contradicted in glo bal simulations performed 
bv lWada fc Normanl (|2007l ) and lWada fc Normanl ([20011 ) 
who do find a lognormal PDF fit without the need for an 
isothermal gas, although note a steepening at high den- 
sities whi ch might be in dicative of a power-law tail on 
the PDF. IKravtsovl (|2003l ) also find a lognormal fit for 
th eir cosmological si mulations of high redshift galaxies as 
do lSlvz et alj (|2005l ) in their 3-D small box simulations, 
but both note that at the high densities where the fit is 
applied, the gas is, in fact, near ly isothermal. 

Bot h the global sim ulations of lWada k, Normanl (|2001[ ) 
and IKravtsovl (|2003h show an insensitivity to the input 
physics, in particular finding that the inclusion of stellar 
feedback does not affe ct the PDF shap e except at the 
low density end where IKravtsovl (|2003h notes feedback 
produces more low density gas. Slyz's small box simula- 
tions at high resolution show a greater sensitivity to the 
inclusion of feedback, without which the gas resides pre- 



dominantly at higher densities. The normalized volume 
weighted gas density PDFs for our simulations are shown 
in Figure [5] at two different times; the left column shows 
the profile at 377 Myrs, when the disk has settled after 
the initial burst of star formation, and the right-hand 
column shows the profile at the end of the simulation, 
when the majority of the gas has been used up in the 
disk. The rows separate our three ISM types with the 
different runs labeled in the key. 

At 377 Myrs, all three of the ISM types show a similar 
shape at high density and are reasonably well matched by 
a lognormal profile. Disks with ISM #1 that do not in- 
clude feedback show a sl ightly bimodal profile, indicativ e 
of a thermal instability (|Vazquez-Semadeni et al.l f2000). 
However, the introduction of feedback smoothes this, al- 
lowing the disk to follow a lognormal profile even at low 
densities. Later evolution of all runs shows an increase in 
the substructure of the disk, with a higher percentage of 
gas found at lower densities as shown in the right-hand 
column of Figure [8] at 1.3 Gyrs. The shape of the PDF 
at this time is similar to that found bv lSlvz et al.l (|2005f ) 
in their small box simulations with feedback, although 
our non-feedback simulations in this ISM show the same 
profile. By this stage, gas depletion is starting to become 
an issue in the runs without feedback whose high den- 
sity tail is being eaten away. This is particularly true in 
the C-type star formation run, C, which shows greater 
gas loss than the D run. Feedback slows this process 
considerably and, although does not totally prevent it, 
it maintains its lognormal fit. The lack of substructure 
and reduced gas loss suggests that feedback is having a 
stabilizing influence on the disk. 

The introduction of background heating in ISM #2 
makes little difference to the profiles at 377 Myrs, pro- 
ducing slightly more low density gas in the feedback 
cases. However, the profile shows markedly less evolu- 
tion than in ISM #1 models, maintaining a larger vol- 
ume at medium densities at the end of the simulation. 
This is most noticeable in the simulations without feed- 
back which have a significantly higher dense gas content 
at t = 1321 Myrs than their ISM #1 counterparts, C and 
D. The amount of low density gas is less than with runs 
with ISM #1, and slightly greater in the feedback case. 
The heating term has the strongest effect on the D-type 
star formation without feedback (HD), reducing gas loss 
by a factor of 10 from the D case so that it resembles the 
feedback runs. This results in the lognormal profile be- 
ing largely maintained in the HCFDBCK, HDFDBCK 
and HD cases. Our higher resolution run, HDHIRES, 
shows little difference from its lower resolution counter- 
part, HD, indicating that these results are not resolution 
dependent. 

The isothermal ISM disks becomes bimodal at around 
700 Myrs, although this is clearly not related to a thermal 
instability. This peak in the low density gas results from 
material that is out of pressure equilibrium in areas such 
the disk center (where star formation has consumed much 
of the gas) and above the disk where we would expect the 
temperature to be high. The log-normal fit is slightly less 
good even at early times here, indicating that a power- 
law fit might b e better. Howeve r , this is very hard to 
judge since like lWada fc Normanl (|2001h . we are limited 
by our resolution. The C- and D-type star formation 
runs show very little difference in this ISM. 
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Fig. 9. — PDF of the volume weighted gas temperature for the simulated galaxy disks over the course of the simulation. 



10* 10° 
T(K) 



10° 



10* 10° 
T(K) 



C 
D 

CFDBCK 
DFDBCK 



HC 
HD 

HCFDBCK 
HDFDBCK 
HDHIRES 



Over the main part of the simulation, our disks show 
a lognormal profile in the medium and high density gas 
which is largely insensitive to ISM model, star formation 
type or the introduction of feedback. This fit stretches 
over several orders of magnitude in density, in keep- 
ing with m odels by I Wada fe Norm"anl d2007l l2001h and 
ISlvz et all (|2005t ) but contradicting |ScaloeialJ ((1991) 
statement that the gas must be isothermal to attain a 
lognormal fit. Substructure in the disk is clearly visi- 
ble in the lower density gas, indicative of the multiphase 
nature of the ISM. 

At late times, the gas profile fit is eroded by gas deple- 
tion in the models where feedback or background heating 
are not present. The addition of background heating pre- 
vents gas depletion for the D-type star formation routine, 
but has less effect on the C-type algorithm. Our pure 
isothermal model shows signs of gas depletion but not as 
marked as in models where cooling is allowed and feed- 
back is not present. The increase of substructure at later 
times makes it difficult to determine the best profile fit. 
The lognormal fit is certainly satisfactory, but a power- 
law curve might well do as good a job, especially in the 
case of the isothermal run which seems to show a steeper 
profile. At our resolution, however, this is still difficult 
to determine. 

4.2. Temperature and Pressure Profile Evolution 

In addition to the volume of gas at different densi- 
ties, the ID temperature distribution can also provide 
valuable information about the ISM. Figure [9] shows the 
evolution of the temperature over the course of the sim- 
ulation. Top panel shows the simulations which have 
ISM #1 while the bottom panel shows ISM #2. The 
isothermal model (ISM #3) has, of course, only gas at 
10 4 K. 

All the simulations with ISM #1 follow a very simi- 



lar evolution, regardless of star formation type or inclu- 
sion of feedback. Initially, the gas quickly cools from 
the 10 4 K starting temperature to the minimum allowed 
temperature of 300 K. As the disk fragments and forms 
stars, the gas is heated by spiral shock waves, infall onto 
the disk and (in the case of runs which include feedback) 
supernovae explosions. By 708 Myrs (third panel), the 
gas is largely in two temperature pockets: warm gas at 
10 4 K and hot gas at 10 6 K. As star formation continues, 
the cooler gas is either used up or heated up, leaving only 
the hotter gas in the disk, which corresponds to the low 
density peak in Figure [5J 

The addition of background heating in ISM #2 re- 
sults in the simulation being significantly more sensi- 
tive to stellar feedback, since the background heating has 
a stronger effect on the no n- feedback runs. The added 
heating stabilizes the warm phase, resulting in the major- 
ity of gas remaining at 10 4 K over the initial part of the 
simulation, rather than cooling to 300 K as in ISM #1. In 
this respect, these disks are most similar to the isother- 
mal case, with the majority of their gas in the stable, 
warm phase. The warm phase remains throughout the 
simulation, but is reduced by the addition of feedback 
which converts more of the gas into the hot phase. This 
is a reflection of what we saw in Figure [5] where HC and 
HD maintained a stronger population of medium den- 
sity gas which we now see is the warm phase, whereas 
the runs HCFDBCK and HDFDBCK resulted in a great 
proportion of gas in the low density (and therefore hot) 
phase. 

Again, the difference between the C- and D-type star 
formation routines is minor, except in the case of HC and 
HD, where HC contains more hot gas that HD. 

The final ID study we can make of this type is the 
volume weighted pressure which is shown in Figure [10] 
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Fig. 10. — PDF of the volume weighted gas pressure for the simu- 
lated galaxy disks after 377 Myrs and end the end of the simulation 
at 1.32 Gyrs. 

for the same times as Figure [SJ While not completely 
isobaric, the range of pressures at 377 Myrs in all the 
non-isothermal disks is reasonably small. This is espe- 
cially true for runs with ISM #2, where radiative heating 
raises the temperature of the lowest pressure gas (we will 
return to this point in more detail in the next section) 
and least true for ISM #3. Feedback acts to broaden 
the distribution. ISMs #1 and #2 show little varia- 
tion over time, implying the disks are in rough pressure 
equilibrium. The isothermal case, however, evolves into 
a bimodal distribution by 1.32 Gyrs, with a significant 
quantity of gas at low pressure. This feature is a direct 
reflection of the density in Figure [5] since the temperature 
cannot change. 

4.3. The ISM as a Three-phase Medium 

Taking t he idea of the tradit i onal t hree-phase ISM por- 
trayed by iMcKee fc Ostriker] (|1977t ) , we next looked at 
the evolution of the cold, warm and hot material in the 
disk over the course of the simulation. Assuming the cold 
ISM consists of temperatures < 10 3 K, the warm ISM is 
between 10 3 K and 10 5 K and the hot ISM is tempera- 
tures above 10 5 K we look at the evolution of both the 
volume and mass fractions in the disks. 

Figure [11] shows the volume weighted (top) and mass 
weighted (bottom) evolution for the disks with ISM #1. 
As in Figure [9] we see the majority of the gas volume 
initially cools, filling the cold phase. As star formation 
removes the dense, cold gas, the warm phase begins to 
dominate. SNe feedback increases the hot phase in the 
CFDBCK and DFDBCK runs, while the non-feedback 
cases become depleted still more of gas, allowing the hot 
phase to dominate which eventually encompasses most 
of the volume in the disk. The resulting pattern is rela- 
tively independent of feedback and star formation type, 
although the move to the warm phase occurs slightly 



earlier for runs with feedback since they have the added 
energy injection. 

The gas mass, on the other hand, resides predomi- 
nantly in the cold phase, sitting in the dense knots of 
matter that go on to form stars. Despite its large vol- 
ume, very little mass actually resides in the hot phase 
(note the change of scale on the hot phase plot abscissa). 

Figure [12] shows the same set of plots for ISM #2. 
The difference between the types of run here is much 
more marked. In all cases however, very little of the 
gas volume is contained in the cold phase, less than that 
for the runs in ISM #1. This is in agreement with Fig- 
ure [5] where we saw a significant proportion of the gas 
stayed at 10 4 K, boosted out of the cold phase by the 
background heating. For the runs without feedback, the 
majority of the volume and mass of the gas sits in the 
warm phase, causing the disk to be closer to isothermal 
in these conditions. This suggests that star formation is 
being suppressed in the HD and HC cases, with respect 
to runs C and D, as the cold star- forming gas is reduced 
by being moved to the warm phase. We will see in sec- 
tion [5] that this is indeed the case. This also causes less 
of the gas volume to be in the hot phase, since the cool 
gas has not all been depleted. The feedback runs also 
start with the bulk of their volume in the warm phase 
but this is swiftly over-ridden by the stronger effects of 
the SNe energy injection, heating the gas and moving it 
into the hot phase. 

A better way to see the phase structure of the ISM is 
via 2D contour plots of the density versus temperature. 
Figure [T2] shows the volume weighted contour plots for 
the disks with ISM #1 (left-hand plots) and disks with 
ISM #2 at 377 Myrs. Straight black lines mark lines of 
constant pressure (neglecting changes in the mean molec- 
ular mass). 

The runs with ISM #1 and no feedback (i.e. C and D) 
show similar behavior: there is rough pressure equilib- 
rium with one or two orders of magnitude scatter along 
the pressure direction (as seen earlier). The sharp fea- 
ture at log(T) = 2.5 is due to our minimum temperature 
cutoff and the gas that resides in the narrow feature that 
rises to higher pressure is gravitationally bound. This 
is the cold, dense phase. A peak in the volume distri- 
bution can be seen for the high-temperature phase (and 
less clearly for the warm phase at T ~ 10 4 K). Feed- 
back acts to broaden this distribution, as can be seen in 
the CFDBCK and DFDBCK runs, making the pressure 
equilibrium less obvious. 

The heated runs (HC and HD) show a very different 
behavior, with a much tighter relationship between den- 
sity and temperature. This is due to the imposition of 
equilibrium between cooling and heating at low tempera- 
tures (below 10 4 K), which leads to the cold phase having 
higher pressure. This tends to force gas into the warm 
phase and reduces the star formation rate because of the 
reduced amount of cold, dense gas. Adding feedback on 
top of this broadens the distribution somewhat but has 
a milder effect than in ISM #1. 

5. OBSERVATIONAL COMPARISON 

So far we have seen that the properties of the ISM 
play a strong role in determining the structure of the 
disk. We have, however, seen some evidence that a uni- 
versal PDF might apply to all disks that could allow the 
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Fig. 11. — Evolution of the cold, warm and hot phases in the disk simulations with ISM #1. Here, the cold phase is defined as gas below 
10 3 K, the warm phase as being between 10 3 K and 10 5 K and the hot phase as having temperatures over 10 s K. The mass fraction refers 
to the gas only. Note different scale for the hot ISM mass fraction. 



Cold ISM: T< 10 3 Warm ISM: 10 3 < T < 10 5 Hot ISM: 1Q 5 < T < 10* 





Fig. 12. — Evolution of the cold, warm and hot phases in the disk simulations with ISM #2. Temperature ranges the same as in Figurc lTll 
Note the different scales on both the cold volume fraction and the hot mass fraction. 



14 



TASKER ANn RRVAN 




?,oo 



3 + 56 
log, J (K) 



3 + 567 
log, J (K) 



Fig. 13. — Two dimensional contour plots for the volume 
weighted gas in the runs marked above. Plots are taken at 
377Myrs. 

observational properties to be largely independent of the 
gas structure. This section focuses on star formation in 
our galaxies and compares the results with the main ob- 
servational relations. 

5.1. Star Formation History 

Initially, we turn to the star formation history to exam- 
ine the rate at which stars are forming over the course of 
our simulations. This is shown in Figure [T4l whose three 
panels depict the evolution of our three ISM types over 
the course of the simulation. In all cases, the curves fol- 
low the same pattern as in TB06; the first 50Myrs sees 
a sudden star burst as the disk becomes gravitationally 
unstable and starts to collapse. The star formation rate 
then peaks and falls off roughly exponentially as gas is 
used up in the disk. 

Considering first the top plot showing disks with 
ISM #1, we see that the runs without feedback have a 
steadily decreasing star formation rate over the majority 
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Fig. 14. — Star formation rate over the course of the simulation 
for the above listed simulations. 



of the simulation. In the last 200Myrs, the simulation 
with D-type star formation starts to show signs of reach- 
ing a constant value, evening out at roughly lM Q yr _1 . 
The C-type run, C, by contrast, heads to zero star for- 
mation and in both these cases, the gas content is signifi- 
cantly depleted as we saw in Figure [51 The time scale for 
this is significantly less than for present day spirals, as 
also found in the lighter disk runs of this type in TB06, 
although of course our initial conditions do not provide 
an accurate cosmological starting point so this is hard to 
judge. The simulations that include feedback, however, 
reach a steady star formation rate at around 600Myrs, 
resulting in a higher star formation rate by 1.3 Gyrs than 
for the non- feedback runs. This shows that feedback 
is causing the star formation to become self-regulated 
which agrees well with the images of the disks in sec- 
tion 13.11 where feedback results the destruction of cold 
clumps of gas. At the end of the run, the feedback sim- 
ulations have a constant star formation rate of around 
10M©yr _1 which is high compared to the Milky Way, 
whose star formation rate is roughly lM Q yr _1 . 

The effect of adding background heating in ISM #2 
causes an overall suppression of the star formation rate. 
In the second panel we can see that the lines are much 
closer together; the non-feedback runs HC and HD hav- 
ing only a slightly lower star formation rate after 1 Gyr 
than the feedback cases. The background heating there- 
fore acts in a similar way to the feedback, suppressing 
the formation of the cold clumps of matter that form 
stars. This is again consistent with Figure O where the 
filamentary structure of the gas was notably less than in 
Figure [TJ Feedback acts to suppress star formation even 
more, with the lines corresponding to HCFDBCK and 
HDFDBCK lying above HC and HD, but the overall dif- 
ference is smaller. Both sets of feedback runs in ISM #1 
and ISM #2 show very similar evolution. At the end of 
the simulation, all the heated runs have a star formation 
rate between 1 - lOM yr _1 . 
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TABLE 2 

Star formation cut-off (radius which includes 99 % of sf) 
and toomre q value at that point. 



Cut-off radius (kpc) Toomre Q Obs. Q 



c 


13.3 


0.40 


0.82 


D 


11.4 


0.26 


0.54 


CFDBCK 


10.7 


0.30 


0.47 


DFDBCK 


12.4 


0.28 


0.58 


HC 


17.8 


6.41 


3.49 


HD 


13.8 


0.95 


0.81 


HCFDBCK 


13.6 


1.02 


0.86 


HDFDBCK 


15.0 


1.45 


1.07 


HDHfRES 


16.1 


1.78 


1.20 


fC 


12.2 


1.83 


0.72 


fD 


12.4 


1.91 


0.76 


ID JEANS 


12.9 


1.87 


0.74 



The two isothermal runs are shown in the bottom panel 
of Figure [HI There is a small distinction between the C- 
and D-type star formation evolution and, contrary to the 
other simulations, the C-type star formation is generally 
higher than the D. This is likely to be due to the C- 
type stars forming in the larger gas clumps than in the 
previous ISM models, due to the increased Jeans length, 
mimicking a behavior closer to the D-type star formation. 
The star formation rates reach a roughly constant value 
after roughly 800Myrs, giving an end rate of lM Q yr -1 , 
close to what we see for the heated cases, HC and HD, 
and in good agreement with the Milky Way. 

5.2. The Star Formation Cut-off 

The gravitational collapse of gas into cold knots of mat- 
ter is the fundamental driving force for star formation in 
our galaxies (this is discussed more thoroughly in TB06 
and references therein). Previous studies, both theoreti- 
cal and observational, suggest that at a given radius the 
density of the gas drops below a given critical value and 
the disk becomes stable, preventing s tar formation from 
occurring beyond this point. Toomre (|Toomrdll964j ) ini- 
tially defined the location of this radius in terms of a sta- 
bility parameter, Q, given by Q = Kc s /-KGY> g , where k 
is the usual epicyclic frequency, c s is the thermal sound 
speed as measured in the disk and S ff is the gas sur- 
face density. Toomre's calculations considered axisym- 
metric perturbations in a single phase two-dimensional 
disk which was found to became gravitationally unsta- 
ble (and therefore able t o for m stars) when Q < 1. 
iGoldreich fc Lvnden-Belll (|1965f ) re-calculated this value 
for a three-dimensional disk to Q < 0.67. Observa- 
tional re sults have also ob served this star formation cut- 
off, with lKennicutt] (|1989h measuring the star formation 
threshold in spiral galaxies to correspond to a Q value of 
1.5. Observations did not allow measurement of the ac- 
tual thermal sound speed of the gas, which was replaced 
by a velocity dispersion of 6 kms -1 . 

In our disks, we measured both the radial cut-off point 
for star formation and the Toomre Q parameter for the 
disk's gas at that radius. We also calculated an observed 
Q value, where we adopt Kennicutt's value of 6kms _1 
instead of the thermal sound speed. It is worth noting, 
however, that these calculations were not entirely sim- 
ple. The greater gas mass and smaller star particle size 
results in a significant gravitational scatter of stars near 
the edge of the stellar disk than we had in TB06, so the 
exact edge of star formation is hard to judge. (This is 



of course true in observational result too, a point we will 
return to in the next section). We therefore take the ra- 
dius where 99 % of the star particles are enclosed as the 
stellar cut-off. The measuring of Q in a multiphase disk 
is an even harder task to perform accurately. Firstly, 
in averaging over an annulus at a given radius, you in- 
clude a wide range of temperatures and densities which 
produce an average Q, not necessary the Q value at the 
star formation sites. Secondly, Toomre's original calcu- 
lations assumed linear perturbations which break down 
at the point of star formation. Ideally, therefore, you 
want to measure the Q value where the disk has become 
gravitationally unstable, but not yet formed stars. In 
TB06, such a point existed, but with our heavier disk 
and smaller star particles the disk begins to form stars 
earlier on. We therefore measure Q at the same time as 
in TB06, approximately 50Myrs after the start of the 
simulation and note that while the disk has fragmented, 
stars have already started to form. 

What is surprising is that we do measure only a small 
range of values for both the cut-off radius and Q as shown 
in Table [2j In all cases, the star formation ends at a 
comparable radius to the Milky Way, which is estimated 
to be around 15 kpc. Within the range we do see some 
patterns. For the non-feedback cases, the C-type star 
formation simulations form stars out to a radius sev- 
eral kiloparsecs greater than in the D-type, something 
we saw visually in Figures [1] and [5] The reason for this is 
that the D-type star formation's higher density thresh- 
old, confining stars to regions where the gas has collapsed 
to form the dense knots of matter. Our C-type star for- 
mation however, can occur at much lower densities, when 
the gas has not fully collapsed. The exception to this is 
the isothermal gas, where, as previously mentioned, the 
larger Jeans length confines the C-type stars to the same 
region at the D-type. The addition of feedback has a 
different effect on the two star formation types. For the 
C-type, it reduces the threshold radius whereas for the 
D-type, the radius for star formation increases. This dif- 
ference is caused by the nature of the resulting outflows 
from the feedback. As we saw in FigureHJ feedback from 
C-type stars produces a smooth, continuous outflow that 
disrupts the dense gas and prevents stars from forming. 
Feedback from the D-type stars however, comes in ener- 
getic bursts which, as we saw in Figure extends the 
gas outflows to much higher radii. These outflows of gas 
can then cool and condense, fragmenting beyond the old 
instability threshold to push star formation out to higher 
radii. 

A slightly surprising result is that the heated disks have 
higher threshold radii than the non-heated cases. Since 
the heating has a stabilizing effect on the disk (as seen in 
Figure fl4|) we would expect it to quench star formation 
in the less dense outer disk. We actually see in Figure Q] 
that the instabilities do extend out further in the DC and 
DD case than the HDC and HDD, but the star forming 
knots stop at a lower radius. This appears to be the re- 
sult of the dissipation of the circular wave that we noted 
occurred in DC in section I5TT1 Although, as mentioned, 
the exact cause of this is unclear. 

If we look at what Q parameter this cut-off corresponds 
to, we find values around 1. Given that Q itself varies 
over several orders of magnitude, the results are very 
uniform, with neither ISM conditions, feedback or star 
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Fig. 15. — Surface star formation rate versus surface gas density 
averaged over the whole disk for different times during the simula- 
tions. This is the Global Sc hmidt law and the solid line shows the 
best fit from observations (Kcnnicutt 1989). 



formation type having a large effect. The only value 
that appears out of place is for run HDC with Q = 3.49. 
Given the range of Q over the disk, this value is still in 
keeping with the other simulations although its higher 
value may be due to problems accurately measuring the 
star formation cut-off. 

Coupled with the lower star formation cut-off, disks 
without heating in ISM #1 have lower Q values than the 
ISM #2 and ISM #3 disks. Feedback and star formation 
type have little effect on these values, which are all below 
1. ISM #2 shows the most sensitivity to stellar condi- 
tions, with the C-type star formation in particular being 
more sensitive to the introduction of feedback than the 
D-type. The isothermal ISM #3 simulations show almost 
no variation with star formation type. For comparison, 
the third column in tableOshows the observationally cal- 
culated Q parameter. The values are reassuringly similar 
to the first calculation of Q, although we can see that for 
ISM #1, the 6kms~ 1 is an overestimate of the thermal 
sound speed, since it raises the value of Q, whereas in 
the other cases it is an underestimate. 

5.3. The Global Star Formation Relation 

Kennicutt's observations of disk galaxies show that the 
surface star formation rate is linked, both on a global and 
local scale, to the surface gas density via Ssfr oc T,^ s . 
In Figure [T51 we show Kennicutt's observational relation 
(solid line) and the results from our different simulations 
for the global version of this law (Global Schmidt law), 
which shows the average surface star formation rate and 
the average gas surface density over the whole disk for 
different (evenly spaced) times over the 1.4Gyrs of the 
simulation. 

Overall, the majority of the simulations reproduce the 
1.5 gradient well, although all over estimate the star for- 



mation rate. The possible exception to this is the isother- 
mal run with D-type star formation, IDD, whose points 
contain a large amount of scatter and a steep decline at 
a gas density of 10M Q pc~ 2 . 

The over-estimate in the SFR was discussed in TB06 
and is likely due, at least in part, to our inability to in- 
clude all the physics at the resolution of the giant molecu- 
lar clouds. Without the added destructive effects of ioniz- 
ing radiation and stellar winds, these clouds lifetimes will 
be prolonged, allowing an unphysically high proportion 
of their gas to be converted into stars. The other possi- 
bility is that neither of our star formation mechanisms is 
very accurate. To a certain extent this is inevitable since 
we d o not resolve the GMC s, but recent work in this 
area (Krumholz & Tan 2006) suggest that while star for- 
mation occurs in the dense molecular clouds, on scales 
larger than a few parsecs the star forma tion efficiency 
is muc h lower than the 50 % suggested by lLada fc Ladal 
( 2003), possible down to 2 %. This reduced efficiency will 
be explored in future papers. 

As in TB06, we find the addition of feedback does act 
to reduce the star formation rate, confirming that we are 
seeing at least the beginning of self-regulation. This is 
most evident in the simulations with C-type star forma- 
tion: C, CFDBCK, HC and HCFDBCK, all of which lie 
close together on the graph. Simulations with D-type 
star formation also follow the observed gradient well, es- 
pecially in the non-heated simulations, D and DFDBCK. 
This is important since, as discussed in TB06, the C-type 
star formation has a Schmidt-like behavior built into it, 
implying that we are simply getting out what we put in. 
However, this is not true for the D-type algorithm, which 
has a simple density cutoff and is otherwise proportional 
to the density. 

As hinted at with the PDFs in Figure [51 the back- 
ground heating in ISM #2 has a stronger effect on the 
HD simulation than on HC, as the higher density regions 
required for this type of star formation are disrupted by 
the increase in temperature. This results in the star for- 
mation at low densities in HD to be significantly reduced 
from D. At these low densities, the results from HD now 
follow the observations very closely, but at higher den- 
sities, the star formation rises above what is expected. 
This suggests that the background heating has less effect 
at the higher densities, where the gas is collapsing despite 
the extra resistance. Introducing feedback reduces this 
effect, still producing an over-estimate of the star forma- 
tion, but the results are now consistent with HCFDBCK, 
CFDBCK and DFDBCK, implying that feedback is the 
dominant effect here. 

The isothermal runs show the poorest agreement with 
observation here. The C-type star formation, IDC, pro- 
duces a gradient slightly steeper than 1.5 for gas densities 
greater than around 6M pc~ 2 but has a large scatter 
below that. 

5.4. The Local Star Formation Relation 

The Schmidt law can also be looked at on a local scale 
where the values are plotted as a function of disk radii 
rather than time. Figure \W\ shows the results for each of 
the simulations where each point represents a different 
radii in the disk. The solid line has a gradient of 1.5, in 
agreement with the observations, but with an arbitrary 
normalization. 
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Fig. 16. — The local Schmidt law shows the variation of surface 
star formation rates with gas surface density at different radii of 
the disk averaged over times between t = 283 - 472 Myrs. The 
solid line is a curve with slope 1.5, as observations indicate (but 
with arbitrary normalization). 

In the top panel where we show the simulations with 
ISM #1, we can see that runs with D-type star forma- 
tion, D and DFDBCK, represent the observed gradient of 
the correlation between the surface gas density and star 
formation rate extremely well. Runs with C-type star 
formation do a reasonable job at high densities, showing 
a larger degree of scatter than the D-type star formation 
runs due to clumps being split by energy injection and 
gravitational interactions and reforming (something that 
is not observed in the D-type run due to the higher bind- 
ing energy of the dense knots). At lower densities, the 
feedback run, CFDBCK, does a noticeable poorer job at 
reproducing the Schmidt law, showing a sharp decrease 
in the star formation rate below gas densities of around 
2M pc -2 , a feature that is not seen in the D-type runs 
which arc not affected by the introduction of feedback. 
This feature is not seen when we introduce background 
heating in ISM #2 (middle panel) where all simulations 
follow the observations very well, although there is again 
noticeably more scatter in the C-type run with feedback 
than in either the run without feedback or the D-type 
simulations, HD and HDFDBCK. 

The isothermal run with D-type star formation follows 
the same pattern, and indeed there is no difference be- 
tween D, HD and ID simulations, showing that the D- 
type star formation is much less influenced by its inter- 
stellar environment. The isothermal run with C-stars, 
however, does not follow the observations well, produc- 
ing a gradient that is much steeper than observed thereby 
underestimating the star formation at low densities and 
over estimating it at high densities. Overall, the range 
in gas densities is much less than in other C-type runs, 
agreeing with the results in section 13.11 that star forma- 
tion is confined to a small section of the disk. 

In both the global and local Schmidt law plots, the 



C-type star formation extends to much lower densities 
than the D-type. Although this is not surprising, since 
its density threshold is lower, it is interesting to com- 
pare this t o observations. A s mentioned in the previ- 
ous section, iKennicuttl () 19891 ) observed a sharp cut-off in 
the star formation in disk galaxies below a critical den- 
sity. More r ecent observations d one using the UV data 
in GALAX (|Boissier et al.l [2006h suggest that this star 
formation cut-off might be a result of the observational 
technique, rather tha n the existence of a critical density. 
iBoissier et ail (|2006l ) looked at 46 spiral galaxies in the 
UV to establish where star formation ended. Previous 
work had been performed by examining the Ha emission 
from galaxies, a technique, this group argues, that makes 
it very difficult to observe low levels of star formation. 
With UV spectra, they found evidence of star formation 
beyond the position of the previously measured cut-off 
radius suggesting that the critical density value is either 
lower than originally measured or possibly does not ac- 
tually exist. If this is the case, then our C-type star 
formation is preferable in this result, since it allows star 
formation to occur at much lower densities in agreement 
with this new result. 

6. DISCUSSION 

What role does the interstellar medium play in deter- 
mining the star formation properties of galaxies? This 
question has two important consequences. The first con- 
cerns the nature of galaxy evolution, asking whether the 
interstellar environment can result in dramatically differ- 
ent structure and star formation properties. The second 
consequence applies to our ability to model galaxy forma- 
tion realistically or whether parsec resolution is needed 
to achieve accurate results. 

In section [3] we visually examined projections of the 
galaxy disks looking at their density, temperature and 
pressure distributions and their star formation. Clear 
differences were seen between our disks, with ISM type, 
star formation algorithm and feedback all playing sig- 
nificant roles. The growth of perturbations in the disk 
was dominated by the conditions in the ISM. In our first 
ISM type, which contained only cooling, we saw the disk 
fragment through the production of a outgoing circular 
wave which collapsed tangentially to form the network of 
filaments and dense knots of gas we saw in Figure [TJ The 
addition of background heating in ISM #2, suppressed 
the formation of the weaker filaments but increased the 
size of the densest knots which extended further out into 
the disk. In our third, isothermal, ISM type, the fila- 
mentary structure was also suppressed and the size of 
the dense knots increased still more, confining their for- 
mation to the central region of the disk. 

These differences made a difference when star forma- 
tion started to occur. In our C-type (cosmological) al- 
gorithm, stars formed over a large fraction of the disk's 
surface, although the largest clusters were confined to 
the densest knots of gas. In the heated case, the heavier 
filaments formed edges around voids of low density, hot 
gas that had occurred through gravitational and thermal 
instability alone. Such voids are reminiscent of porous 
ISMs seen in the HI maps of some dwarf galaxies. In 
ISM #1, further fragmentation had prevented this from 
occurring. The isothermal disk confined the stars to the 
central disk region where the large knots formed and the 
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resulting heavier clusters scattered the smaller ones out 
of the disk. The extent of this scattering was extreme, 
and care should be taken if an isothermal model is used 
with this star formation recipe to avoid unphysical re- 
sults. The effect in ISM #3 of confining the star for- 
mation to the central region meant that there was little 
difference between the C-type algorithm and the D-type 
star formation algorithm where we restricted the forma- 
tion of stars to the densest clouds. This was investigated 
quantitatively in section [5.21 where we measured the star 
formation cut-off radius and compared it to the Toomrc 
stability criterion Q at that point. In the ISM #3 case, 
the cut-off radius was almost identical in both the C- and 
D-type star formation routines whereas in ISM #1 and 
ISM #2, the C-type stars extended further into the disk. 
The Q stability criteria itself was found to be around 1, in 
good agreement with both analytical and observational 
results. 

This changed again when we included feedback. In 
both the ISM #1 and ISM #2 cases, feedback destroyed 
the dense knots of gas, suppressing the star formation as 
we saw in Figure 1141 Its addition largely wipes out the 
structural differences between the two ISM types and 
their profiles, both in the projections and the vertical 
profiles in Figures [6] and [7] follow similar shapes. Feed- 
back also causes gas to be ejected from the disk's surface 
in a galactic fountain, something the isothermal disk is 
incapable of replicating. The different star formation al- 
gorithms have a much larger effect on the feedback than 
ISM type. The D-type algorithm, confining the star for- 
mation to the densest knots, focuses the energy injection 
from the feedback, causing the outflows to be bursts in 
different areas of the disk. The C-type, by contrast, acts 
equally across the disk and at a more continuous level. 
This has two effects on the disk. The first was seen in 
Figure 0] where the outflows from the disk's surface are 
uneven. The second is seen in the star formation cut-off, 
where the feedback in the D-type case triggers fragmenta- 
tion of the gas in the stable regions of the disk, extending 
the star formation threshold. 

Section 2] further examines the ISM by looking at the 
evolution of its density and temperature over time. We 
see that background heating stabilizes a significant pro- 
portion of the gas in the warm phase in Figure QTJ bring- 
ing it nearer to the isothermal state as can be seen in the 
2D contour plots in Figure [131 This could explain the 
success of models using an isothermal ISM; it is a gross 
simplification of the real system, but if a single phase be- 
comes stabilized in the disk, an isothermal gas can be a 
good approximation. The addition of feedback, however, 
dramatically increases the gas volume in the hot phase 
causing the heated and non-heated disks to show a wide 
variety of continuous phases. Both ISM #1 and ISM #2 
produce disks that are largely isobaric, in keeping with 
analytical expectations. The addition of feedback acts as 
a pressure value to eject gas from the disk, ensuring the 
disk remains isobaric across its height. The isothermal 
disk, by contrast, has strong pressure variability since the 
fixed temperature means that the pressure has to reflect 
the density distribution. 

In section[5]we turned to look at the observable proper- 
ties of our disks, including star formation history, cut-off 
and the relation between gas surface density and surface 
star formation rate (Schmidt laws). No one model pro- 
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Fig. 17. — Two dimensional contour plots for the volume 
weighted gas in the run with increased heating source. Plots are 
taken at 377 Myrs. 



duces these results significantly better than any other. 
A result that implies, despite structural differences, that 
global star formation is not strongly dependent on the 
interstellar environment. The disks overall produce the 
observed properties reasonably well, having a star forma- 
tion threshold of between 11-18 kpc in agreement with es- 
timations for the Milky Way and producing the Schmidt 
observed gradient of 1.5. The slight exception to this 
were the isothermal cases which reproduced the observa- 
tional result poorly at low densities. Of the two star 
formation algorithms, each has its own strengths and 
weaknesses. The C-type algorithm allows the extension 
of the Schmidt law to much lower densities, where recent 
observational evidence suggests there is star formation. 
However, the small stellar clusters gravitationally inter- 
act to be scattered in the disk, making it more difficult 
to determine the cut-off density for star formation and 
increasing the scatter in the local Schmidt relation. The 
D-type algorithm, meanwhile, show a tight correlation 
in the local Schmidt relation but appears to reproduce 
the gradient of the global relation less well. However, it 
is also the least affected by ISM type, since it consid- 
ers only the densest gas structures. Possibly a hybrid 
of these two schemes would combine their strengths to 
produce the best results. 

As mentioned in section 12.11 we performed an addi- 
tional run with a higher background heating source. This 
run was identical to run HD, but with the background 
heating increased by a factor of 30. The main difference 
the additional heating had on the disk was to increase 
the pressure of the ISM. This occurred because the en- 
hanced heating rate moved the gas in thermal equilib- 
rium to higher temperature at fixed density, pushing the 
contours shown in Figure [T3l to the right and up. This 
can be seen in Figure [17] The effect this has on the 
star formation, however, was fairly minimal. There was 
a slightly increased star formation rate at early times 
which led to a greater gas depletion near the end of the 
simulation. This resulted in the global Schmidt law fol- 
lowing a relation identical to that seen in the HD run, 
but at a slightly raised star formation rate. 

Ultimately, HI images of the whirlpool galaxy and 
LMC tell us that the ISM plays an important role in 
the galaxy's evolution and it must be born in mind that 
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while we are able to resolve a multiphase structure in 
our ISM, we are no yet at the detail of the small box 
simulations. However, the uniformity of our results for 
the star formation history and Schmidt laws suggest that 
the exact details of the gas may be simplified and still 
achieve a correct star formation production. 

7. CONCLUSIONS 

We performed high resolution simulations of global 
disk galaxies with three different interstellar medium 
properties: our first model allowed radiative cooling of 
gas, the second model allowed radiative cooling and in- 
cluded a background heating term while the third model 
held the gas at a constant temperature. For each of these 
ISM models, we considered two different star formation 
mechanisms; a 'C type that allowed low density forma- 
tion of stars at a correspondingly low efficiency and a 
'D' type that only allowed star formation in the densest 
clumps, but with a high efficiency. We also investigated 
the addition of feedback from type II supernovae for non- 
isothermal disks. Overall we found: 

1. The structure of the disk is strongly affected by the 
ISM environment. Background heating suppresses 
small-scale perturbations and increases the size of 
the star forming knots. This additionally helps to 
regulate star formation and reduce gas depletion 
in the disk. The isothermal equation of state in- 
creases the Jeans' length in the disk and leads to 
the formation of much larger clumps (probably un- 
physically large) which results in extreme gravita- 
tional scattering of nearby star particles. It also 
confines star formation to the central region, re- 
gardless of the critical density specified in the star 
formation routine. This leads to gas depletion in 
the disk center. 

2. The addition of feedback destroys star- forming 
clumps and causes gas to be ejected off the disk's 
surface. The nature of the outflows is dependent on 
the star formation algorithm with the low threshold 
density, C-type routine producing a more continu- 
ous flow compared to the D-type strong interme- 
diate bursts. Both cases result in star formation 
being suppressed in the disk as dense knots of gas 
are destroyed, but the feedback in the D-type case 
results in triggered star formation in the outer re- 
gions of the disk. 

3. The structure of the ISM in both ISM #1 and 
ISM #2 show a continuous range of densities and 
temperatures that are not well represented by a 
single phase model. The addition of heating re- 
duces the range of values found in the non-feedback 
case, in particular stabilizing the warm phase at the 
same temperature as the isothermal ISM. However, 
the introduction of feedback significantly increases 
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the range of densities and temperatures, almost 
eliminating the signature of background heating. 

4. The PDFs for all disks are well represented by 
a lognormal curve over several orders of magni- 
tude. This is largely insensitive to the introduc- 
tion of stellar feedback or ISM environment. At 
later times (and for the isothermal run at earlier 
times) a power-law fit would also be possible and 
it is hard to tell, at current resolutions, which of 
these two fits would work best. 

5. The star formation in the disks are shown to cut- 
off at a radius comparable to the Milky Way when 
the Toomre Q parameter is around 1. Exact mea- 
surements are difficult to achieve in a multiphase 
medium. 

6. All simulations reproduce the slope of the observed 
relation between star formation and gas surface 
density well on both the global and local scale. The 
possible exception is the isothermal run with C- 
type star formation which shows a steeper drop off 
in star formation rate on local scales. There is some 
evidence that D-type star formation does not pro- 
duce the gradient as well at C-type on global scales 
and that its low density cut-off may be higher than 
recent observations. 

Overall, we conclude that the interstellar medium is a 
multiphase environment strongly affected by background 
heating, feedback and cooling and plays a significant 
part in the disk's structural evolution. It appears, 
however, from these preliminary calculations that its 
exact structure can be simplified and still achieve the 
correct star formation properties for disk galaxies. The 
exception to this is possibly the isothermal disk which 
did not produce a realistic multi-phase model, lead to 
very large clump formation, and can never produce a 
galactic fountain effect. 



EJT and GLB would like to thank Adrianne Slyz 
for helpful comments and suggestions and the referee 
for helping to improve the presentation of the pa- 
per. EJT acknowledges support from a Theoretical 
Astrophysics Postdoctoral Fellowship from Dept. of 
Astronomy/ CL AS, University of Florida. Both au- 
thors acknowledge the National Center for Supercom- 
puting Applications and the University of Florida High- 
Pcrformance Computing Center for providing computa- 
tional resources and support and also support from NSF 
grants AST-05-07161, AST-05-47823, and AST-06-06959 
that have contributed to the research results reported 
within this paper. 



Bryan,G.L. Comp. Phys. and Eng. 1999, 1:2, p. 

B ryan, G. L., No rman, M. L. & Abel, T. 2001, arXiv: 

astro-ph/0112089 proceedings of Supcrcomputing 2001 
Cen, R. & Ostriker, J. P. 1992, ApJ, 399, L113 
Cox, D. P. 2005, ARA&A, 43, 337 
Deul, E. R., & den Hartog, R. H. 1990, A&A, 229, 362 



20 



TASKER AND BRYAN 



Elmegreen, B. G. 2002, ApJ, 577, 206 

Ferguson, A. M. N., Irwin, M. J., Ibata, R. A., Lewis, G. F., & 

Tanvir, N. R. 2002, AJ, 124, 1452 
Gill, S. P., Bryan, G. L., & Tasker, E. J. 2007, in prep. 
Goldreich, P. & Lynden-Bell, D. 1965, MNRAS, 130, 97 
Harfst, S., Theis, C, & Hensler, G. 2006, A&A, 449, 509 
Joung, M. K. R., fe Mac Low, M. -M. 2005, ArXiv Astrophysics 

e-prints, arXiv:astro-ph/0601005 
Kennicutt, R. U. 1989, ApJ, 344, 685 
Kennicutt, R. C. 1998, ApJ, 498, 541 

Kim, S., Staveley-Smith, L., Dopita, M. A., Freeman, K. C, Sault, 

R. J., Kesteven, M. J., & McConnell, D. 1998, ApJ, 503, 674 
Kravtsov, A. V. 2003, ApJ, 590, LI 

K rumholz, M. R., &: Tan, J . C. 2006, ArXiv Astrophysics e-prints, 

|arXiv:ast ro-ph/0606277 
Lada, (J. J., h Lada, E. A. 2003, ARA&A, 41, 57 
Li, Y., Mac Low, M-M., & Klessen, R. 2005, ApJ, 620, 19L 
Li, Y., Mac Low, M.-M., & Klessen, R. S. 2005, ApJ, 626, 823 
MacLow, M.-M. 2004, Ap&SS, 289, 323 
Marri, S. & White, S. D. M. 2003, MNRAS, 345, 561 
McKee, C. F., & Ostriker, J. P. 1977, ApJ, 218, 148 
Navarro, J., Frenk, C, White, S.D.M. 1997, ApJ, 
Norman, M. L. & Bryan, G. L. 1999, ASSL Vol. 240: Numerical 

Astrophysics, 19 
Norman, C. A., & Ikeuchi, S. 1989, ApJ, 345, 372 
O'Shea, B. W., Bryan, G., Bordner, J., Norman, M. L., Abel, T., 

Harkness, R., & Kritsuk, A. 2004, ArXiv Astrophysics e-prints, 

|astro-ph/ 0403044 
Otte, B., Dixon, W. V., & Sankrit, R. 2003, ApJ, 586, L53 
Robertson, B., Yoshida, N., Springel, V., & Hernquist, L. 2004, 

ApJ, 606, 32 



Rhode, K. L., Salzer, J. J., Westpfahl, D. J., & Radice, L. A. 1999, 

AJ, 118, 323 
Rosen, A. & Bregman, J. N. 1995, ApJ, 440, 634 
Sarazin, C. L. & White, R. E. 1987, ApJ, 320, 32 
Scalo, J., Vazquez-Semadeni, E., Chappell, D., & Passot, T. 1998, 

ApJ, 504, 835 
Schmidt, M. 1959, ApJ, 129, 243 
Semelin, B. & Combes, F. 2002, A&A, 388, 826 
Slyz, A. D., Devriendt, J. E. G., Bryan, G., & Silk, J. 2005, 

MNRAS, 356, 737 
Springel, V. & Hernquist, L. 2003, MNRAS, 339, 289 
Stone, J. M. & Norman, M. L. 1992, ApJS, 80, 753 
Tasker, E. J., & Bryan, G. L. 2006, ApJ, 641, 878 
Tasker, E. J., & Bryan, G. L. 2006, ApJ, 642, L5 
Toomre, A. 1964, ApJ, 139, 1217 

Truelove, J. K., Klein, R. I., McKee, C. F., Holliman, J. H., Howell, 

L. H., & Greenough, J. A. 1997, ApJ, 489, L179 
van der Hulst, J. M. 1996, ASP Conf. Ser. 106: The Minnesota 

Lectures on Extragalactic Neutral Hydrogen, 106, 47 
Vazquez-Semadeni, E., Gazol, A., & Scalo, J. 2000, ApJ, 540, 271 
Wada, K., Spaans, M., & Kim, S. 2000, ApJ, 540, 797 
Wada, K. & Norman, C. A. 2001, ApJ, 547, 172 
Wada, K., Meurer, G., & Norman, C. A. 2002, ApJ, 577, 197 
Wada, K., & Norman, C. A. 2007, ApJ, 660, 276 
Wolfire, M. G., Hollenbach, D., McKee, C. F., Tielens, A. G. G. M., 

& Bakes, E. L. O. 1995, ApJ, 443, 152 
Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, 

A. G. G. M. 2003, ApJ, 587, 278 



